Small- World Bonds and Patchy Percolation on the Hanoi Network 
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The bond-percolation properties of the Hanoi networks are analyzed with the renormalization 
group. Unlike scale-free networks, they are meant to provide an analytically tractable interpolation 
between finite dimensional, lattice-based models and their mean-field limits. In percolation, the hi- 
erarchical small-world bonds in the Hanoi networks impose a new form of order by uniting otherwise 
disconnected, local clusters. This "patchy" order results in merely a finite probability to obtain a 
spanning cluster for certain ranges of the bond probability, unlike the usual 0-1 transition found on 
ordinary lattices. The various networks studied here exhibit a range of phase behaviors, depending 
on the prevalence of those long-range bonds. Fixed points in general exhibit non- universal behavior. 
PACS: 64.60.ae, 64.60.ah, 64.60.aq 
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INTRODUCTION 



Percolation — the formation of a large connected 
component — is a geometric property in the arrangement 
of a many-body system that can strongly impact its phys- 
ical behavior [l|, 0]. Typically, to observe any form of 
emergent, collective phenomena, an extensive fraction 
of the degrees of freedom must be sufficiently intercon- 
nected. Examples abound, from transport in amorphous 
materials [3| or the onset of plasticity in complex fluids [4] 
to the spreading of rumors or disease in social networks. 
In particular, while percolation on regular and fractal 
lattice geometries [H and ordinary random graphs 0, 0] 
is a well-developed industry, its properties and impact of 
transport on the many conceivable forms of engineered 
networks (random, scale-free, etc) are just beginning to 
be explored 0, [a, @|. 

Here, we study the percolation properties of the re- 
cently introduced set of Hanoi networks (if! [TT1 . IT3 |. 
These networks mimic the behavior of small world sys- 
tems without the usual disorder inherent in the construc- 
tion of such networks. Instead, they attain these proper- 
ties in a recursive, hierarchical manner that lends itself to 
exact renormalization [T]|. These networks do not pos- 
sess a scale-free degree distribution [3,[ll,[TB|; they are, 
like the original Watts-Strogatz "small worlds" [17], fig , 
of regular degree or have an exponential degree distribu- 
tion. Yet, in turn the Hanoi networks have a more "phys- 
ically" desirable geometry. In particular, they have the 
potential to provide an analytically tractable interpola- 
tion between finite dimensional, lattice-based models and 
their mean-field limits (although no such interpolation 
will be considered here [HI].) Scale-free networks geared 
towards social phenomena or other lattices often used to 
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analyze physical systems, such as the hierarchical lattice 
derived from the Migdal-Kadanoff bond-moving scheme 
figl . [20I . [2l| , do not have a physically relevant mean-field 
behavior. 

Many of the phenomena of interest in statistical 
physics, such as critical behavior, are fundamentally re- 
lated to their percolation properties. We find that the 
renormalization group (RG) applied to the Hanoi net- 
works behaves very differently from the traditional ap- 
proach but similar to a low-dimensional hierarchical lat- 
tice with small-world bonds that was recently introduced 
[H, [HI • This is surprising, as this behavior has usually 
been associated with scale-free graphs. Due to the hierar- 
chical nature of the small- world bonds, the RG possesses 
unrenormalized parameters, which change the charac- 
ter of the RG-flow in unusual ways. In the exactly ob- 
tained phase diagrams, stable and unstable fixed points 
are drawn out on lines that often merge in branch points, 
generally making the scaling near fixed points parameter- 
dependent and, hence, non-universal. As the flow can get 
attracted onto a stable line, entirely new, mixed ordered- 
disordered phases appear in what we call "patchy" or- 
der: Otherwise isolated patches of localized clusters at 
low levels of the hierarchy get reconnected globally at 
higher levels with finite probability (depending on said 
parameters) to attain long-range order. In the percola- 
tion problems discussed here, this manifests itself in the 
fact that for a fixed bond-probability the infinite network 
may percolate with a finite probability between zero and 
unity. (In conventional percolation, that probability has 
an instant 0-1 transition at the critical bond density p c , 
although the size of the spanning cluster varies continu- 
ously, see Fig. 17 in Ref. [H). 

Our paper is structured as follows: In the next sec- 
tion, we review the construction of the Hanoi networks, 
followed by a pedagogical discussion of the renormaliza- 
tion group applied to percolation on hierarchical lattices 
with small- world bonds in Sec. IIIII Then, in Sec. |IV] we 
study the percolation properties of the Hanoi networks 



Figure 1: Depiction of the 3-regular network HN3 on a semi- 
infinite line. The entire graph becomes 3-regular with a self- 
loop at n — 0. Note that HN3 is planar. 

in some depth, and we conclude in Sec. [Vl 



Figure 2: Depiction of the planar network HN5, a comprised 
of HN3 (black lines) with the addition of further long-range 
bonds (green-shaded lines). There is no distinction made be- 
tween black and shaded lines in our studies here. 



II. GEOMETRY OF THE HANOI NETWORKS 

Each of the Hanoi network possesses a simple ge- 
ometric backbone, a one-dimensional line of sites 
n, either infinite (— 2 k < n < 2 fe , k — ► oo), semi-infinite 
(0 < n < 2 k , k — > oo), or closed into a ring of N = 2 k 
sites. Most importantly, all sites are connected to their 
nearest neighbors, ensuring the existence of the ld- 
backbone. To generate the small-world hierarchy in these 
networks, consider parameterizing any integer n (except 
for zero) uniquely in terms of two integers via 

n = ?(2j + l). (1) 

Here, i > denotes the level in the hierarchy whereas 
j labels consecutive sites within each hierarchy. For in- 
stance, i — refers to all odd integers, i = 1 to all in- 
tegers once divisible by 2 (i. e., 2, 6, 10,...), and so on. 
In these networks, aside from the backbone, each site is 
also connected with (one or both) of its nearest neighbors 
within the hierarchy. For example, we obtain a 3-regular 
network HN3 (best done on a semi-infinite line) by con- 
necting first the backbone, then 1 to 3, 5 to 7, 9 to 11, 
etc, for i = 0, next 2 to 6, 10 to 14, etc, for i = 1, and 4 
to 12, 20 to 28, etc, for i = 2, and so on, as depicted in 
Fig. [TJ (The corresponding 4-regular network HN4 [Icj ]. 
which we will not study here, is obtained in the same 
manner but connecting to both nearest neighbors in the 
hierarchy, i. e., 1 to 3, 3 to 5, 5 to 7, etc, for i = 0, 2 to 
6, 6 to 10, etc, for i = 1, and so forth.) 

Previously [lfl], it was found that the average chemical 
path between sites on HN3 scales with the distance I as 

d Hm ~ Vl (2) 

along the backbone. In some ways, this property is rem- 
iniscent of a square-lattice consisting of N lattice sites, 
whose diameter (=diagonal) is also ~ \/~N. 



While HN3 (and HN4 P3) are of a fixed, finite degree, 
we introduce here convenient generalizations of HN3 that 
lead to new, revealing insights into small-world phenom- 
ena. First, we can extend HN3 in the following manner to 
obtain a new planar network of average degree 5, hence 
called HN5, at the price of a distribution in the degrees 
that is exponentially falling. In addition to the bonds in 
HN3, in HN5 we also connect all even sites to both of 
their nearest neighboring sites that are within the same 
level of the hierarchy i > 1 in Eq. ([1]). The resulting 
network remains planar but now sites have a hierarchy- 
dependent degree, as shown in Fig. [21 To obtain the 
average degree, we observe that \ of all sites have de- 
gree 3, \ have degree 5, | have degree 7, and so on, 
leading to an exponentially falling degree distribution of 
V {a = 2i + 3} oc 2~\ Then, the total number of bonds 
L in the system of size N = 2 k is 

fc-i 

2L = 2k - 1+ ( 2i + 3 ) 2 ^ 1_i = 5 2 fc - 6, (3) 

8=0 

and thus, the average degree is 

9 T 

(a) = — ~ 5. (4) 

In HN5, the end-to-end distance is trivially 1, see 
Fig. [21 Therefore, we define as the diameter the largest of 
the shortest paths possible between any two sites, which 
are typically odd-index sites furthest away from long- 
distance bonds. For the N = 32 site network depicted in 
Fig. [21 for instance, that diameter is 5, measured between 
site 3 and 19 (starting with n — is the left-most site), 
although there are many other such pairs. It is easy to 
show recursively that this diameter grows strictly as 

d HN5 = 2[k/2\+l~log 2 N. (5) 
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Figure 3: Depiction of the non-planar Hanoi network HN-NP. 
Again, starting from a ld-backbone (black lines), a set of long- 
range bonds (blue-shaded lines) is added that break planarity 
but maintain the hierarchical pattern set out in Eq. JTJ. The 
RG on this network remains exact. 

We have checked numerically that the average shortest 
path between any two sites appears to increase logarith- 
mically with system size N as well. 

The networks HN3 and HN5 have the convenient but 
(from a mean-field perspective) unrealistic restriction of 
being planar. In fact, with a minor extension of the def- 
inition, it is easy to also design Hanoi networks that are 
both non-planar and fully renormalizable. The simplest 
such network, which we dub HN-NP here, is depicted in 
Fig. El (Extensions will be introduced in Ref. [24j].) 

To obtain the average degree, we observe that \ + 4 
of all sites have degree 3, | has degree 5, ^ has degree 
7, and so on, leading to an exponentially falling degree 
distribution, as for HN5. The total number of bonds L 
in the system of size iV = 2 k is 

fc-l 

2L = 3 2 k - 1 + J2{2i-l)2 k -' 1 + 2k+(2k-2) 

i=2 

= 4 2 fe -4, (6) 
and thus, the average degree is 

= 4. (7) 

Here, too, it is easy to see that the shortest paths between 
sites increases logarithmically with system size N . 

III. RG FOR PERCOLATION 

The principles of the renormalization group (RG) for 
percolation as it will be use in this paper are easiest to 
explain for a simple hierarchical lattice [2l|, |22| . (Percola- 
tion on these hierarchical lattices has recently also been 
discussed in Ref. [13] •) I n this recursively constructed 
lattice, bonds from a preceding generation g are built up 
in the next generation g + 1 according to a fixed pattern, 
until the complete network is constructed at a generation 




Figure 4: Pattern for the recursive generation of a hierarchical 
lattice. In zeroth generation (top), the lattice consists of a 
single bond between the two "boundary" sites (squares). In 
successive generations each bond of that type (solid line) is 
replaced by a new type of bond (dashed line) in parallel with a 
sequence of two of the original bonds (solid lines), see middle 
and bottom for generation one and two, resp. 

g = k — > oo. Applying real-space renormalization to the 
complete network in effect reverses the recursive built-up 
of the network. The complete network with the "bare" 
operators obtained in the g — fc-th generation provides 
the initial conditions, or the n = 0-th step in the RG, 
which we evolve to n — > oo (i. e., k — > oo) and study the 
fixed points of this recursion and their stability. 

An example of this is shown in Fig. [4] for the first three 
generations (g = 0, 1, 2) of a hierarchical lattice that cor- 
responds to a one-dimensional line (solid bonds) with a 
hierarchy of small- world bonds (dashed lines) [23|. It 
is fruitful to regard both types of bonds as distinct at 
this stage, although in any complete network one may be 
interested mostly in the case where all bonds are indistin- 
guishable and exist with equal probability in the context 
of percolation, say. 

In Fig. we list all 8 graph-lets consisting of those 
three bonds. Summing the weights of those 5 graph- 
lets that connect end-to-end and, thus, contribute to the 
percolation probability q n +i, we obtain 

q n+1 = p + (1 -p)ql- (8) 

It should be noted that there is no corresponding recur- 
sion for the small- world bond p (dashed line). It enters 
the RG anew at each step, bare and unrenormalized. (As 
we will see, this feature is a defining characteristic also 
of RG in the Hanoi networks, as discussed in Ref. [10].) 

It is very simple, or course, to obtain the fixed points 
of the RG in Eq. |(8]) in the thermodynamic limit k — > oo 
with the Ansatz q n +i — q n — q* for n — > oo. Our case 
yields the trivial unstable fixed point q* = 1 and a non- 
trivial but stable fixed point at 

**Cp) = t^ (p<\)- < 9 > 

We have obtained a line of stable fixed points depending 
on the parameter p. For p = 0, we obtain the ordinary 
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Figure 5: All possible combinations of reducing bonds from 
the n-th RG step to form a new bond at step n+ 1, reversing 
the pattern set in Fig. [4] All graphs on the left percolate (i. e., 
connect end-to-end) by merit of the dashed bond alone, while 
on the right only the top graph provides such a connection. 
If q„ is the probability for a solid and p for a dashed bond, 
the weights for the graphs on the left are pq^, pq n (1 — q n ), 
pq„ (1 — q„), andp?n (1 — q-nf" and on the right are (1 — p)q'i, 
(1 -p)q n (1 - q n ), (1 -p)q n (1 - q n ), and (1 - p)q n (1 - q n ) 2 
from top to bottom. 



and trivial percolation problem of a one-dimensional lat- 
tice with a percolation threshold at q c — 1: For any initial 
bond probability qo < 1, the RG flow from the recursion 
in Eq. J5]) evolves away from the unstable fixed point 
q* = 1 to the stable fixed point q* — 0. The ordered 
(percolating) state is attained only for qo = 1. For all 
p > 0, the phase diagram for the RG flow is shown in 
Fig. [6j For any p < h , the usual ordered state remains 




0.1 0.2 0.3 0.4 0.5 0.6 



P 

Figure 6: Phase diagram for the RG flow of the simple hi- 
erarchichal lattice due to Eq. l|8]). For any given p, the flow 
for q n evolves along the flow trajectories starting from a cho- 
sen initial value qo- Except for points (p < \,qo — l) on the 
(thick red) line of unstable fixed points, all points (p, qo < 1) 
flow towards the (thick black) line of stable fixed points from 
Eq. J9} or q* = 1 for p > ^. This is especially true for a ho- 
mogeneous choice of bond probabilities in a network, qo = p 
along the (dash-dotted) diagonal line. 



confined to the unstable fixed point qo = q* = 1. For any 
other starting value < qo < 1, the RG flow converges 
onto the line of non-trivial stable fixed points given by 
Eq. |(9|) in which the usual non-percolating, disordered 
state is now replace by a partially ordered state, i. e., 
even in the thermodynamic limit there is a finite prob- 
ability to percolate. For p>\, the stable and unstable 
fixed points merge and an ordered, percolating state is 
reached irrespective of qo (even when it vanishes) , merely 
by the strength of the small-world bond, which repeat- 
edly, at any generation of the hierarchical lattice, has a 
chance p to rectify any non-percolating sub-lattices from 
previous generations (see, for instance, the bottom net- 
work in Fig. 21). In reference to the corresponding effect 
we found previously for the Ising ferromagnet on HN5, 
we call this partially ordered state "patchy" [2J|: ordering 
occurs between the lucky neighborhoods that happen to 
be connected by the occasional small-world bond. 

We obtain this picture from a local analysis of Eq. JH]) 
near the fixed point with the Ansatz 

q n ~ q* - 5 n (10) 

assuming 5 n <C 1- To leading order with q* = 1 yields 

6 n+1 ~ 2(1 -p)6 n (11) 

with solution 

S n ~ [2{l-p)] n S . (12) 

For 2 (1 — p) < 1, the corrections 8 n contract for increas- 
ing n, leaving the fixed point stable. In contrast, for 
2 (1 — p) > 1 the corrections escalate (eventually violat- 
ing the assumption) and the fixed point is said to be 
unstable. As shown in Fig. [6J the line of unstable fixed 
points at q* = 1 for small p thus disappears at p — \ and 
is replaced by a stable fixed point at all larger p. 

Typically, in a complete network, we make no dis- 
tinction between the probabilities of any of the bonds 
in the network, whether they are small- world or not. 
Hence, in light of Fig. El the RG flow would initiate 
with qo = p, which corresponds to the flow trajectories 
emanating from the diagonal (dash-dotted) line. Con- 
sequently, projecting the flow in the diagram just onto 
that diagonal, one would conclude that there is no con- 
ventional phase transition. Starting on the open diagonal 
at any point the flow proceeds to a larger effective per- 
colation probability. There is an interesting transition in 
this behavior at p — \ between a fully ordered (perco- 
lating) phase above, where this flow always converges to 
q* = 1, and a partially-ordered, patchy phase below, con- 
verging to some non-trivial value of q* = q*(p = qo) < 1 
given by Eq. ([9]) . We will see that this situation resembles 
very closely that of percolation on HN5. Notice that the 
unstable fixed point discussed above, typically itself the 
focus of any study in critical phenomena, seems to have 
become irrelevant to these considerations, as our initial 
conditions never cross an unstable manifold. 
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Figure 7: Pattern for the recursive generation of a 2d hier- 
archical lattice, originating from the Migdal-Kadanoff bond- 
moving scheme in a square lattice, with the addition of a 
small-world bond. Each single, solid bond between the two 
sites (squares) of a previous generation is replaced by a small- 
world bond (dashed line) in parallel with two sequences of two 
of the original bonds (solid lines) at the next generation. 



It is, in fact, just as simple to analyze the correspond- 
ing behavior for the small-world hierarchical network on 
an underlying two-dimensional lattice [12] . As shown in 
Fig. [3 we merely need to add a second sequence of two 
(solid line) bonds to the previous hierarchical generation 
of the network in Fig. [U resulting in a new graph-let of 
five bonds. As in Fig. \E[ we consider the now 2 5 = 32 
graph-lets regarding the end-to-end connectivity and sum 
the weights of those that contribute to percolation. Anal- 
ogous to the one-dimensional RG recursion ([8]) we find 



In+x = p + (1 - (2 - t£) 



(13) 



which is very similar in character: it also depends on a 
non-renormalized, bare parameter p for the probabilities 
of the small-world bond at every recursion. The same 
fixed-point analysis as above leads to a more complex but 
straightforward algebraic equation for q* = q*{p), and 
we plot the respective RG flow in Fig. [8l Without small- 
world bonds, at p = 0, we already observe the emer- 
gence of a nontrivial unstable fixed of the ordinary two- 
dimensional hierarchical lattice at q* = l/4> = 0.618 . . ., 
where <f> = (\/5 + l) /2 is the "golden ratio" [H- In 
this case, q*(p) develops a remarkable branch-point sin- 
gularity within the physical domain at the critical point 
P = -§2, where two branches of fixed points pinch off into 
the complex plane and the stable fixed point attained for 
all points along the diagonal qo = p jumps discontinu- 
ously from q* (p — > p~) — | to q* (p — > p + ) = 1. 

Again, we can verify these findings with a local analysis 
along the line of fixed points with the same Ansatz as in 
Eq. ifTOl , where q* (p) is now a solution of Eq. (jT3j) , see 
the parabolic line in Fig. [H We obtain 



4(l-p)q* l-(q 



(14) 



where the simultaneous solution of the fixed point equa- 
tion together with the stability condition on the Ansatz, 




Figure 8: Phase diagram for the RG flow of the 2d hier- 
archichal lattice due to Eq. (| 1 3 j) . using the same notation 
as in Fig. [6] Here, there exists an additional phase for 
p < p — j|j = 0.15625 due to the (thick red) line of un- 
stable fixed points on some line qo(p) < 1- Points above 
that line flow into the now stable ordered fixed point q* = 1. 
Note that for the homogeneous choice of bond probabilities 
in a network, i. e., go = p along the (dash-dotted) diagonal 
line, we observe a discontinuous behavior at £>: Coming from 
the partially ordered phase below, the flow towards the finite 
percolation probability q* < 1 jumps (along the thick dashed 
line) towards q* = 1 in the fully ordered phase above p. 



defines the branch of stable fixed points, while its vio- 
lation implies the unstable branch. Both branches meet 
when the inequality in Eq. lfT5|) is saturated, at p = ^ 
and q* = |. The fixed point analysis for q* — 1 leads to 



S n 



+i 



4(1 -p) 5$ 



This non-linear relation has the solution 



4(l-p) 



[4(1 -p) 5 ] 



(16) 



(17) 



4(l-p)q* l-(q*) 



< 1, 



(15) 



Deriving from our assumption, it is certainly So < \ and 
hence S n vanishes and the fixed point proves to be exceed- 
ingly attractive and stable, ever more so for increasing p. 



IV. RG FOR PERCOLATION ON HANOI 
NETWORKS 

In this section, we study percolation on HN3, HN5 and 
HN-NP with the renormalization group. It is easy to con- 
vince oneself that percolation on HN3 can only occur at 
full connectivity p — > 1, because HN3 is finitely ramified: 
To prevent end-to-end connectivity or the emergence of 
a giant component only a fixed number of bonds need to 
be cut at any system size [H]. This is most apparent by 
plotting HN3 as a branched Koch curve, see Fig. [9l Yet, 
the limit p — > 1 in HN3 itself is interesting, showing log- 
arithmic finite-size corrections, which appear impossible 
to resolve numerically. Moreover, the derivation of the 
RG recursions for HN3 are almost identical to that for 
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Figure 9: Representation of HN3 as a branched Koch curve. 
The one-dimensional backbone shown in Fig. [l] is marked in 
black here; long-range bonds are shaded in red. In this repre- 
sentation, the diameter in Eq. ([2j corresponds to the base-line 
in the Koch curve. To prevent end-to-end connectivity, only 
two bonds right down the symmetry axis need to be cut, in- 
dependent of size. Hence, the network is finitely ramified. 



A 



B 



C 



HN5 and HN-NP, which in turn have a more interesting 
phase diagram. 



A. Renormalization Group Studies for HN3 

First, we consider the renormalization group for end- 
to-end percolation on HN3. In this case, all steps can 
be done exactly. Following Taitelbaum et al. [26j] in 
their analysis of the Sierpinski gasket, the RG step con- 
sists of eliminating every second (odd) site in the net- 
work and calculating new, renormalized probabilities for 
connections between the remaining sites, as explained in 
Fig. HOl Four independent probabilities (R n , S n ,T n ,U n ) 
are evolved under RG; the fifth probability, N n , follows 
from conservation of probability. 

We obtain a system of RG recursion equa- 
tions, in which the renormalized probabili- 
ties (i?„+i, S n +i, T„ + i, U n +i) are functions of 
(R n , S n ,T n ,U n ;p). As for the case of the hierar- 
chical lattices in Sec. IIII1 these equations are unusual, as 
they retain a memory of the original bond probability p 
for all times, which enters as a source-term th roug h the 
dependence on the next level in the hierarchy flOj. The 
initial conditions can be read off the right graph-let in 
Fig. [TUl Unlike in the other sections, where we merely 
execute a local stability analysis near fixed points, here 
we need to make contact with the initial conditions in 
the analysis to explore finite-size corrections. Therefore, 
we follow again the notation in Ref. [26] and shift 
the RG recursion to evolve from initially n = — k to 
terminate at n = for a system of size L — 2 k in 
the thermodynamic limit k — > oo. Originally there is 
no bond ac, thus U-k = and R-k — p 2 because it 
only consists of having both, ab and be, connected. We 
have S-k = T-k = p(l — p) for an exclusive ab or be 
connection, respectively, hence N-k = (1 — p) 2 ■ 

We determine the form of the recursion equations for 
these probabilities by a simple counting procedure, as 



Figure 10: Depiction of the (exact) renormalization group 
step for percolation on HN3. The step consists of elimi- 
nating every second site (with label b) in the upper plot 
and expressing the renormalized connection probabilities 
(Rn+i, S n +i, T n +i, U n +i, N n +i) on the lower plot in terms of 
the old values (R n , S„,T n ,U„, N n ). Here, R n refers to the 
probability that three consecutive points abc are connected, 
S„ refers to ab being connected but not c, T n to be being 
connected but not a, U„ to ac being connected but not b, 
whereas N„ refers to the probability that neither a, b, or c 
are mutually connected. Upper case letters A,B,C refer to 
the renormalized sites at n + 1. Note that the original net- 
work in Fig.[T]does not contain bonds of type U„ (long-dashed 
line), but that they become relevant during the RG process. 



Rn+l 
S n +i 



explained in Fig. [TTJ We find 

R 2 n + 2U n R n + U 2 + 2 P S n R n , 

+U n S n + pS n T n + (1 -p)R n S n , 
R n T n + R n N n + U n N n + U n T n 

+U n S n +pS n T„ + (1 -p)R n S n , 
pS 2 , 



T, 



n+l 



U n +1 

N n+1 



(18) 



T 2 + 2N n S n + 2N n T n + N 2 

a 



+2(1 -p)S n T n + (1 -p)S 2 n . 

Note that unlike for the Sierpinski gasket [26] but similar 
to Eqs. ((8]) and (fl3|l . the RG recursions depend explic- 
itly on p at every step, not merely through the initial 
conditions. 

Obviously, the equations imply that T n = S n for all n, 
making the T-operator redundant. Combining probabil- 
ities S n and T n into a single probability S n for having 
a connection between either ab or be, and rescaling by a 
factor of h, Eqs. (fl8ll can be simplified to 



Rn+l — (Rn + U n ) + pS n Rn, 
S n +1 = 2 (S n + U n ) (Rn + U n ) 

Un+l = ^rpS 2 , 



n^n pRnS n , 



(19) 
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for all n. The initial conditions at n 
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Figure 11: Example of graph- lets contributing to the RG step 
in HN3. Eliminating every second site in the network (dots la- 
beled b), the connectedness of the renormalized sites {A, B, C) 
have to be assessed. There are 7 bonds in the elementary 
graph-let: 2 each for pair-connections ab, be, and ac, and 
one previously unrenormalized bond with raw probability p 
stretching between the two sites b. Bonds of higher level in the 
hierarchy (vertical lines) remain unaffected until a later RG 
step. Thus, 2 7 = 128 graph-lets with combinations of these 
bond must be considered. Only four examples are displayed 
here. In the first, all bonds are present with a direct bond 
between abc on both sides of B; each such triangle has a prob- 
ability -Rn/4, and the bond p is present, making the weight 
pRn/16. Because, after renormalization, the remaining sites 
ABC are also directly linked, this graph-let contributes to 
Rn+i- In the second example, only ac but not 6 are linked 
on the left, contributing a weight of U n , but there is a path 
connecting abc on the right with weight -Rn/4, and there is no 
bond p, making the total weight (1 —p)U„R n . Again, there is 
a path between ABC, so it contributes to Rn+i- In the third 
example, the S n bond between ab on the left and the T„ bond 
on the right together with the p-bond ensure a connection be- 
tween AB, thus making a contribution to S„+i. The 4th 
example is the only graph-let that contributes to U n +i with 
a direct link between AC but not B. 



N n+ i = (S n + N n ) 2 -~pSl, 



which satisfies the constraint 



R-k 
S—k 

U- k 



P , 

2p(l-p), 
0, 



(21) 



Note that the fixed point equations can be solved ex- 
actly for arbitrary p, as they involve the solution of three 
coupled quadratic equations, but the expressions are ex- 
tremely complex. Instead, we search for fixed points of 
these equations numerically. For any physical p inside 
the unit interval we find no physical fixed points with 
< N*,R*,S*,U* < 1. 

The analysis of the RG recursion Eqs. I|19p proceeds 
as follows. The obligatory trivial fixed points, i. e., 
stationary solutions of Eqs. US]), at (N* , R* , S* ,U*) = 
(1,0, 0, 0) and (0, 1, 0, 0) exist for any p. We search for a 
non-trivial fixed point by choosing a value of < p < 1 , 
and numerically evolve Eqs. lfT9|) for that value of p, 
starting from the initial conditions in ([2~T]l . We find that 
for all values of p, the equations evolve towards the stable 
fixed point at (1,0,0,0), characteristic for a completely 
disconnected network, although this evolution takes ever 
longer, the closer p is to unity. Therefore, we conclude: 



Pc 



= 1. 



(22) 



This result is not too surprising, because this network 
is of finite ramification: a bounded number of cuts can 
separate the network into two extensive pieces. Finitely 
ramified objects, such as the triangular Sierpinski gas- 
ket, can only percolate with certainty at full connectiv- 
ity, p c = 1. Corrections to this behavior in HN3 are more 
interesting and, in fact, decay faster for k — » oo than for 
other finitely ramified objects like the Sierpinski gasket. 



1. First- Order Correction to p = 1: 

To determine the finite-size corrections for p c — 1, we 
expand around the unstable fixed point (0, 1, 0, 0) to first 
order with the Ansatz 



Rn 



s n l \ (■si 1) «l), 



„(i) 

'In i 



(•tf J « 1 



(23) 



Inserting into Eqs. I|19p and dropping terms beyond first 
order, we obtain immediately that 



= 4P = o. 

The norm in Eq. lj20j) implies 1 ~ 1 — r„ + Sn \ i. e., 



1 = N n + R n + S n + U n 



(20) 



.(i) 



(24) 
(25) 
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and we get from the equations in (fT9l) for S n +f. 



3 (D 

'n+l 



(26) 



which together with Eq. f25|) have the consistent solution 



,(1) _ o(l) 



A(2- P y\ 



(27) 



Applying the initial conditions in f2"Tj) , we can now assess 
the finite-size corrections to p c = 1 with the Ansatz 



1 



£fe, 



(fife « 1) ■ 



Inserting this and the result for and s^jj, in Eq. 
into Eqs. ([21)1 yields a consistent result, which is 



(28) 

G3 



A(l + e fc ) A 



2e* 



(29) 



Note also that with p in Eq. (|28|) . the last initial condition 
in Eq. J27J) evaluates consistently to N-k = (e 2 ) . 
To leading order, we obtain from Eq. lj29j) 



In A; 



(30) 



independent of the unknown constant A. In general it 
would be best to extract e& directly from Eq. ((29j) for 
any given k, as Eq. lj30j) is poorly convergent. Unfortu- 
nately, though, we have not found a way to fix A in this 
local expansion. We conclude here by noting that the 
leading corrections to p c — 1 as stated in Eq. f30|) decays 
faster than what Ref. [26j found for the Sierpinski gasket, 

ej; ~ 1/ ^2-\/A;J , where k — log 2 iV. We continue with the 
second-order correction in the Appendix, with the final 
result given in Eq. (|57|l . 



B. RG for HN5 

Just as HN5 is a simple extension of HN3, see Sec. IH 
above, so is the analysis of the percolation problem. Yet, 
the added bonds in HN5 precisely rectify the finite ramifi- 
cation that made the percolation phase diagram for HN3 
trivial. On the face of it, it is not entirely obvious that 
percolation on HN5 would be any more interesting, be- 
cause its ramification is only marginally infinite. Now 
the number of lines to be cut grows logarithmically with 
size N — 2 k , because each new level k of the hierarchy 
adds another end-to-end connection, see Fig. [2j 

The RG step for HN5 involves merely one extra bond 
in the elementary graph- let, see Fig. [T2l The deriva- 
tion of the RG recursion equations is similar to Eq. (fl9| , 
however, we now have to evaluate twice as many graph- 
lets. But one half of the graph-lets — those without the 
new line — simply correspond exactly to those of HN3 in 
weight (times a factor of 1 — p) and in the operator to 
which they contribute. The other half again have the 
same weight (times a factor p, for the now-present new 
bond), but the additional bond changes the operator to 



P 
P 



by c 



^1 



\b 
— O— 



\ 



A 



B 



A 



B 



C 



Figure 12: Depiction of the (exact) RG step for percolation 
on HN5. In contrast to the corresponding step for HN3 in 
Fig. 1101 the long-dashed bonds are germane to the network, 
not emergent properties. Thus, the long-dashed bond on top 
has still probability p; the long-dashed bonds of the previous 
hierarchical level below have a probability of p plus the emer- 
gent contribution from the previous step. The graph-let after 
the RG step (bottom) is identical to that in Fig. 1101 



which these graph-lets contribute in 69 out of the 128 
cases by adding extra connectivity. These considera- 
tions and the corresponding discussion along the lines 
of Sec. [IVA] lead to 

Rn+l = (Rn + U n f + 2p (S n + N n ) (R„ + U n ) 

p 2 

+p(l-p)S n Rn + y5 n , 



S n +i — 



2(l-p) [(S n + N n ) (R n + U n ) 

-p(l -p) R n S ni 



U, 



n+l 



N, 



n+l 



(Sn + N n 



(31) 



which again satisfies the norm-constraint in Eq. ([20"ll . 
The initial conditions are now 



Ro = p 2 (3 - 2p) , 
So = 2p(l-p) 2 , 

U 
N 



p(\ -py 
(i-p) 3 - 



(32) 



The algebraic fixed-point equations obtained from 
Eqs. I|3ip for n«B|l-»co are extremely complex. It 
is easy to establish the fixed point at (N* , R* , S* ,U*) — 
(0,1,0,0) for any p. But the other trivial stable fixed 
point at (N*,R*,S*,U*) = (1,0,0,0), describing the dis- 
ordered state, only exists at p = 0. Instead, we find (nu- 
merically, as the solution of a sixth-order polynomial) a 



p 



Figure 13: Phase diagram for the RG flow of HN5 due to 
Eq. (131 [I . The all-connected operator R is used as an indicator 
of percolative order in the system. For any given p, R n evolves 
along the shown trajectories starting from an initial value 
Ro- Except for points (p < p = 2 — <j>, Ro = 1) on the (thick 
red) line of unstable fixed points, all points (p, Ro < 1) flow 
towards the (thick black) line of stable fixed points obtained 
from Eq. l[3~T]l or to R* = 1 for p > p. For a homogeneous 
choice of bond probabilities in the network given in Eqs. l132[| . 
i. e., Ro = p 2 (3 — 2p) along the dash-dotted line, the stable 
fixed point is always approached from below. 

nontrivial equation for a line of fixed-point which merges 
into the fully ordered, percolating fixed-point R* = 1 
that is stable for all p > p — 2 — cj> = 0.381966 . . ., where 
<t> = (v5 + l) /2. The resulting phase diagram in Fig. [13] 
is very similar to that in Fig. [6] for the hierarchical small- 
world lattice. The difference is that HN5 is more mean- 
field like in the sense of ordinary random graphs, with an 
exponential instead of a scale-free degree distribution. 

We obtain the value for p explicitly from the analysis 
of the fixed point at R* = 1. In fact, we proceed similar 
to Sec. IIV Al by applying the same Ansatz in Eq. (|38j) to 
the RG recursions in J37|) . A first-order expansion yields 



r (i) = s (i) 
'n+l 



(2 -p) 



(33) 



with Un ' = T]n^ = 0. Finally, we have 



_ a (i) 



A[(2-p)(l-p)Y 



Despite the obvious similarity to Eqs. (|38H27p . the key 
difference of the solution in Eq. ((34j) is that there is a non- 
trivial transition at p — 2 — (f>, derived from the marginal 
stability condition, 



1 



(2-p)(l-p) 



(35) 



such that the fixed point is unstable for < p < p 
and stable for p < p < 1. Clearly, the approach to 
the fixed point is highly parameter-dependent, leading 
to non-universality such as in the correlation-length ex- 
ponent for p < p: [l3| 

" = ta[(2-pKi-p)r (36) 



A 



B 



C 



Figure 14: Depiction of the (exact) RG step for percolation 
on HN-NP. In contrast to HN3, the long-range bond connect- 
ing bb in Fig. [10] is replace by a crossing pair of such bonds 
connecting ab. The £/-bond connecting ac is an emergent op- 
erator, like in HN3. The graph-let after the RG step (bottom) 
is identical to that for HN3 in Fig. \M and HN5 in Fig. [12] 



C. Percolation in HN-NP 

A more substantial extension to the previous planar 
Hanoi networks is given by HN-NP introduced in Sec. H] 
see Fig. [3] It extends the basic network-building concept 
into the non-planar regime while preserving exact renor- 
malizability. The single, long-range bond of HN3 in the 
elementary graph-let in Fig.[l0]is here replaced by a pair 
of crossing small- world bonds, as explained in Fig.HH As 
the second diagram shows, the graph- let resulting from 
the RG step is identical to those for HN3 and HN5. As 
in HN3, only the [/-operator that connect sites ac ex- 
clusively emerges anew during the RG recursion. (Else- 
where, we also consider that operator as pre-existing for 
HN-NP, the same step that converts HN3 into HN5 pi].) 

Evaluation of the 2 s = 256 elementary graph-lets as 
described in connection with Fig. Qj] above results in the 
recursion equations 



(34) Rn+i = (Rn + U n y +p(3-p)R n S n + 2pR n N n 



-pS n U n + -p s. 



2 c2 



(2 - p) [(1 - p) R n + U n ] S n + 2 (1 - p) R n N n 
+2U n N n + P N n S n + P (l-p)S 2 n , 



P ( 1 ~ \p ) Sn+pN^S,. 



U n +i 

N n+1 = Nl + 2{l-p)S n N n + {l-pfS 2 n , 



(37) 



obeying the normalization constraint in Eq. (|20]1 . and 
with the same initial conditions as for HN3 in Eqs. (|2l)l . 

At first glance, the phase diagram for HN-NP appears 
to have "cured" all the unusual small- world effects of the 
previous networks and reverted to ordinary mean-field- 
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Figure 15: Phase diagram for the RG flow of HN-NP due to 
Eq. H37II . Here, the all-connected operator R is used as in 
indicator of percolative order in the system. For any given p, 
R n evolves along the shown trajectories starting from an ini- 
tial value Ro- The homogeneous choice of bond probabilities 
in the network given in Eqs. IpTjl . i. e., Ro = p 2 along the 
(dash-dotted) diagonal line, crosses the line of unstable fixed 
points just below the branch-point singularity. 



like behavior: There is a domain for small p > where 
the percolation probability vanishes and one at larger 
p < 1 where percolation is certain, with what seems to 
be the usual 0-1 transition in the probability at some crit- 
ical p c (see, e. g., Fig. 17 in Ref. |l|). Yet, closer inspec- 
tion of a relatively narrow regime near that transition, 
about for 0.3 < p < 0.4, reveals a more subtle behavior. 
If we didn't have the exact RG equations in (|37|) avail- 
able, a numerical simulation may well have missed it due 
to finite-size corrections. Within that narrow region, en- 
larged in Fig.mi the resulting phase diagram is similar to 
Fig. O with a branch-point where two lines of stable and 
unstable fixed points merge and disappear into the com- 
plex domain. But this branch-point is oriented in the op- 
posite direction, with a (minute) jump in the percolation 
probability at pi = 0.319445181 . . . preceding a continu- 
ous rise towards, again, p u — 2 — tfi = 0.381966 . . . above 
which percolation is certain. Unlike for Fig. [U R* = 
and R* = 1 are always fixed points (variably stable or 
unstable) for all p: While in the 2d hierarchical lattice 
the small-world bond for any p can boot-strap even the 
empty network into a spanning cluster with finite proba- 
bility, in HN-NP this is only possible for a disconnected 
network when p > p m = | . 

The analysis near the fixed points provides more details 
about these observations. The results for the fixed point 
at R* = 1 for p < p u is identical to that for HN5 in 
Sec. IIV B[ see Eq. l(33|) . and provides again p u = 2 — 4> 
as the transition point between the line of unstable and 
stable fixed points. In turn, the analysis near R* = 0, 
is novel. Enforcing the norm-constraint N n = 1 — R n — 
S n — U n and expanding around the unstable fixed point 
to first order with the Ansatz 



R,, 



< 1 



57, 

Un 



«<*) 



(tf « l) , 



(38) 



leads to 



'n+l 
'n+l 



"n+l 



(1) 



2pr, 

2(l-p)rW+psW+2u£\ 



ps 



(i) 



The first relation immediately yields 
r« ~ (2p)"r« 



(39) 



(40) 



and the physically relevant solution for the remaining 
coupled set of linear recursions provides 



.4 



P 



1 



1 



P 



(41) 



(The modulus of the second, oscillatory, solution decays 
with n for all < p < 1.) Again, the highly non-trivial 
scaling exponent derived from Eq. l[4"T]) is non-universal. 
For small p < p m = | , all corrections decay and the fixed 

point remains stable. For p > p m , the solutions for Sn 
and Un" 1 first become divergent, enough to make the fixed 
point itself unstable. Yet, initially remains conver- 
gent, slowing the RG flow. For p > \, the fixed point 
rapidly becomes unstable. The analysis of the branch- 
point at pi can only be done numerically because it in- 
volves the solution of sufficiently high-order polynomials 
to defy an explicit treatment. 



V. CONCLUSIONS 

We have obtained the exact phase diagrams for bond- 
percolation on variants of the Hanoi networks using RG. 
For a small fraction of long-range bonds and a regular 
degree distribution in HN3, which leaves the network 
finitely ramified and with an average distance that grows 
like a square-root of the system size, percolation requires 
a full bond-density, p c = 1. The hierarchical structure 
there manifests itself only in the finite-size corrections. 
Adding more long-range bonds at the price of obtain- 
ing an exponential degree distribution but a mean-field 
like logarithmic distance between sites results in the net- 
work HN5. There, it leads to the unusual phenomenon 
of percolation with a finite probability, which we dubbed 
"patchy" order, facilitated by the addition of ever more 
long-range bonds at higher levels of the hierarchy that 
can reconnect patches of localized clusters at lower levels 
of the hierarchy. Introducing a non-planar version of a 
Hanoi graph provides nearly ideal mean-field behavior, 
with a completely disordered phase at low bond-density 
and a fully ordered phase at high density. Yet, even in 
this case a narrow window with mixed behavior inter- 
venes. 
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Here, we have all but scratched the surface of the anal- 
ysis of percolation on small- world networks. In light of 
the ever increasing importance of small-world networks, 
practical and conceptual, we will focus on exploring fur- 
ther characteristics, such as details of clustering behav- 
iors and correlations, on the variety of network configu- 
rations represented here. Only then can we hope to find 
solid classifications on the possible behaviors found in 
the real world. For instance, it might be possible to find 
recently observed discontinuous transitions in the size of 
spanning clusters ("explosive" percolation [2?], [H, HH E3] 
also in the neighborhood of the branch points of these 
merely hierarchically constructed networks. 



Appendix 

Second-Order Correction to p c = 1 in HN3: 



The inhomogeneous term in this relation can be neglected 
at this order, in fact, as the expression in the square 
brackets vanishes in the limit p — ► 1, and we are left with 
the same solution for sl 2 ^ as in first order, 



a (2) 



B(2-p) n , 



(46) 



where B<i, (In fact, we expect B/A ~ e^.) 

f 2") 

Regarding the relation for r„ in Eq. (f43]) , we can drop 
the last term in this order as 1 — p — > and find 



„( 2 ) 

r n+l 



2r (2) +2u ( 2 ) +ps ( 2 ) ; 

2rl 2 > + P -A 2 (2- p) 2n - 2 + pB{2 — p) n .(47) 



The homogeneous solution of this relation at n = — k 
would decay exponentially, ~ 2~ fc , and can be discarded, 
leaving us with 



We expand around the unstable fixed point (0, 1, 0, 0) 
of the RG recursions (fT9ll for HN3 to second order with 
the Ansatz 



R n ~ 1-rW+ri 2 ), (rf « r« « l) , 
Sn ~ rW+ S L 2) , (4 2) «r««l), 
U n ~u£\ (4?«r^«l), (42) 
N n - V^, (^«rW«l), 

where we have used Eq. lj25|) to eliminate s„ in favor of 

■ Of course, all terms of lower order cancel, and after 
dropping terms of higher than second order, we arrive at 



(2) 
r n+l 

J2) 
'n+1 



(2) 
U n+l 

(2) 
Vn+1 



2r^+2u^+p S ^ + (l-p) (r,«) 2 , 
(2-p) S ^+2^+^-p-2^ (V«) 2 , 

5 (*>)'' (43) 



l--p)(r. 



(i) 



assuming that squares of terms of first order are of second 
order. The solutions for and rjn^ are 

ui? = \A\2- v r-\ 

A 2 (2-p) 2 - 2 . (44) 

(2) 

Inserting these results into the equation for s\L\ yields 



,(2) 



(2) 
S n+1 



(2~P)4 2) 



l-(2- P ) 2 A 2 {2-p) 



(45) 

2n-2 



i 2) ~ fA 2 (2- p) 2 - 4 



- P B(2-p) n+1 . (48) 



Inserting the second-order solutions in Eqs. (|44|46I48[) 
together with the results from Sec. IIV ATI into Eqs. f42|) 
and evaluating at g = k, we get 



R-k 
S-k 

U- k -i 
N, 



l-A(2-p) k +pB(2-p) 
+P-A^2-py 2k -\ 

A{2- P y k + B{2- P y\ 
\A\2-py 2 \ 

X-\p\ AH2~py 2k -\ 



-k-l 



(A > B), 
(49) 



We extend the correction on p c to second order, i. e., 

p ~ l-e k + be 2 kl (50) 

where b = b(ek) should only vary weakly with e k . In- 
serted into Eqs. lf49|) . we get 



R-k 



l-2e k + 4 



2e fe 



S-k 

U-k+i 
N- k 

where we used 



-26 In 



2 + 2&ln 
2e fe 



2efe 
A 



B 
> 



A 



B 
> 

'ekA 



-fe> 



; -ki 



(51) 



A(2- P y 



- A(l + e k -be 2 k ) , 
~ A(l + e k y k [1 + bkel 

~ 2e k 1 - 6 In (jjj^) e k 
because Eq. ((29j) yields -ke k ~ hx{2e k /A). 



(52) 



12 



Correspondingly, we can expand the initial conditions 



R-k = P 2 

~ l-2e fe + e£[l + 26] 

5_ fe = 2p(l-p) 

~ 2e k + 4 [-2 - 26] , 

7V_ fc = (1-p) 2 



(53) 



Note that in Eq. (|5T|) we have used the result for the 
newly emerging operator U at g = k — 1, because U- k 
itself would strictly vanish to all orders. But with the 
help of its relation in Eq. lfl9|) , we obtain via the initial 
condition for S n : 



fc+i 



= P 3 a-pY 



(54) 



4 



Comparison of those results with Eqs. I|5ip prove the 
consistency for N and U. As expected from Sec. IIV A 1[ 
the relations for R and S are also consistent immediately 
up to first order. Requiring consistency to second order 
provides relations 



26 
2b 



2 + 26 In 



2tu 
A 



} B 



-26 In 



2eu 
A 



t B 

! ^4' 



that determine the integration constant B and, more im- 
portantly, the second-order correction 6 to p c : 



6 = -i 



1 



(55) 



B = — e k A. 
4 

A redundant test of these results is provided by the nor- 
malization constraint, 



1 = R. k + S- k + M fc , 



(56) 



which we find satisfied to second order. (Remember that 

U- k = 0.) 

Hence, we obtain finally for the finite-size corrections 
of p c : 



1 

1 - £fc — 7 



1 



4 1 + K^) 



4 



(57) 



where the relation in Eq. lj29| . A(l + e k )~ ~ 2e k , de- 
fines efc. Unfortunately, we do not see any way to de- 
termine the arbitrary constant A in this expansion. We 
assume that it is a reflection of the local nature of our 
expansion for k — > oo [HI]. It seems that A would have 
to be determined from global information about the so- 
lution, such as from an evaluation at n = 0. 
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